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We investigate the temperature-disorder (T-a) phase diagram of a three-dimensional gauge glass 
model, which is a cubic-lattice nearest-neighbor XY model with quenched random phase shifts Axy 
at the bonds, by numerical Monte Carlo simulations. We consider the uncorrelated phase-shift 
distribution P{Axy) ~ exp[{cosAxy) / a], which has the pure XY model and the uniform distribution 
of random shifts as extreme cases at cr = and cr ^ oo respectively, and which gives rise to equal 
magnetic and overlap correlation functions when T = a. 

While the high-temperature phase is always paramagnetic, at low temperatures there is a ferro- 
magnetic phase for weak disorder (small a) and a glassy phase at large disorder (large a). These three 
phases are separated by transition lines with different magnetic and glassy critical behaviors. The 
disorder induced by the random shifts turns out to be irrelevant at the paramagnetic-ferromagnetic 
transition line, where the critical behavior belongs to the 3D XY universality class of pure systems; 
disorder gives only rise to very slowly decaying scaling corrections. The glassy critical behavior 
along the finite-temperature paramagnetic-glassy transition line belongs to the gauge-glass univer- 
sality class, with a quite large exponent f — 3.2(4). These transition lines meet at a multicritical 
point M, located at Tm = o"a/ = 0.7840(2). The low-temperature ferromagnetic and glassy phases 
are separated by a third transition line, from M down to the T = axis, which is slightly reentrant. 

PACS numbers: 75.50.Lk,05.70.Fh,64.60.F-,74.62.En 



I. INTRODUCTION AND SUMMARY 

Spin glass models are simplified, although still quite 
complex, models retaining the main features of physical 
systems which show glassy phases. They may be con- 
sidered as theoretical laboratories where the combined 
effects of quenched disorder and frustration can be stud- 
ied. Their phase diagrams and critical behaviors provide 
examples of possible scenarios which can be used to inter- 
pret the experimental results of complex materials. For 
example, Ising-type spin glasses, such as the ±J Ising 
model model disordered uniaxial magnetic materi- 
als characterized by random ferromagnetic and antiferro- 
magnetic short-ranged interactions. While many theoret- 
ical and numerical works have been devoted to the study 
of the phase diagrams and magnetic and glassy critical 
behaviors of Ising- like spin glasses, see, e.g., Refs. [2-11 and 
references therein, much less is known about the thermo- 
dynamic properties of spin glass models with continuous 
symmetries. In the case of the Heisenberg spin glass there 
is some evidence for a finite-temperature glassy transi- 
tion, but its nature is still debated, see, e.g., Refs. H 
and references therein. 

Another physically interesting spin glass model is the 
XY model with random shifts, also known as the gauge- 
glass model, which is characterized by a global U(l) sym- 



metry. It has been proposed as a simplified model of 
disordered granular superconductors, to describe vortex- 
glass transitions (THlSj. The phase diagram and criti- 
cal behaviors of two-dimensional (2D) gauge glasses have 
been much investigated, see, e.g., Refs. [ij, and ref- 
erences therein. It is now well established that at weak 
disorder there is a low-temperature quasi-long-range or- 
der phase separated by a Kosterlitz-Thouless transition 
line from the paramagnetic phase; no glassy phase ex- 
ists at finite temperature, but a T = glassy criti- 
cal behavior at sufficiently large disorder. A discus- 
sion of the general temperature-disorder phase diagram 
of the three-dimensional (3D) gauge glass can be found 
in Refs. [H, [O- The ferromagnetic phase of the pure 
cubic-lattice XY model is expected to survive at weak 
disorder while it should disappear at large disorder, 
where a low-temperature glassy phase may exist. Sev- 
eral numerical works have addressed the existence of a 
vortex-glass phase at finite temperature in the case of a 
uniform random-shift distribution p^j - lsij . providing ev- 
idence of a finite-temperature glassy transition. Apart 
from these results for the extreme disordered case, the 
phase diagram and the critical behaviors at the different 
transition lines have not been numerically investigated 
yet. We also mention that experimental results for vor- 
tex glassphases in superconductors have been reported 
in Refs. 133 - 1361 but the experimental scenario for the be- 
havior at the transition does not appear settled yet. 
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In this paper we investigate the temperature-disorder 
phase diagram of the 3D gauge glass and its critical be- 
haviors along the transition lines which separate the dif- 
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FIG. 1: Sketch of the temperature-disorder {T-a) phase dia- 
gram of the 3D gauge glass. 



FIG. 2: The dashed hnes sketch the T and a values of our 
MC simulations. 



ferent phases. 

The 3D gauge-glass model is defined by the partition 
function 



Z{{A}) = J [Dyj] expi-n/T), 



(1) 



■H = - ^ Reip^Uxyipy = - ^ cos(6'^ - (^y - ^xy), 

(xy) (xy) 

where T is the temperature, ip^ = e*^^, Uxy = 6*"^"=", 
and the sum runs over the bonds {xy) of a cubic lat- 
tice. The phases A^y are uncorrelated quenched ran- 
dom variables with zero average. A Gaussian distribu- 
tion Pa{Axy) oc exp(— A^j^/(2cr)) is often considered in 
the studies of gauge glasses. In this paper we consider a 
slightly different cosine distribution 



P{Axy) OC exp 



/ coaAj^y 



(2) 



Analogously to the Gaussian distribution, we recover the 
pure cubic-lattice XY model for cr and uniformly 
distributed random phase shifts in the limit cr — > cx). The 
cosine distribution is particularly interesting because it 
lends itself to some exact relations along the so-called 
Nishimori (N) line T = a [l^ , such as the equality 
of the magnetic and overlap correlation functions, which 
are useful to identify the multicritical point where the 
transition lines meet in the T-a phase diagram. 

Our study of the T-a phase diagram, and the mag- 
netic and glassy critical behaviors along its transition 
lines, is mostly based on numerical Monte Carlo (MC) 
simulations. Supplementing the numerical results with 
renormalization-group (RG) and finite-size scaling (FSS) 
analyses, we arrive at the phase diagram sketched in 
Fig. [TJ In Fig. [2] we show where we performed the MC 
simulations in the T-a plane. Our main results are the 
following. 

While the high-temperature phase is always paramag- 
netic, at low temperatures we have a ferromagnetic phase 



for weak disorder (small a) and a glassy phase at suffi- 
ciently large disorder (large a). These three phases are 
separated by different transition lines: a paramagnetic- 
ferromagnetic (PF) transition line, a paramagnetic- 
glassy (PG) transition line and a ferromagnetic-glassy 
(FG) transition line, meeting at a multicritical point M 
located along the N line. The phase diagram of 3D gauge 
glasses presents several analogies with the temperature- 
disorder phase diagram of 3D ± J Ising spin glasses 
where analogous phases and transition lines appear. 

The 3D gauge glass shows different magnetic and 
glassy critical behaviors at the transition lines separating 
the different phases. 

We argue that the disorder induced by the random 
shifts is irrelevant at the PF transition line starting from 
the pure XY transition point for cr = 0, at [131 Txy = 
2.201842(5). Thus the asymptotic critical behavior be- 
longs to the 3D XY universality class of pure systems 
with U(l) symmetry, characterized by the correlation- 
length critical exponent 37] i^xy — 0.6717(1). Analo- 
gously to randomly dilute 3D XY models , the disor- 
der induces new scaling corrections which gets suppressed 
very slowly, as 0(^~'^'* ) where f is the critical length scale 
and cjd = 3 - 2/vxY = 0.0225(5). A FSS analysis of 
MC simulations at cr = 0.35, up to L = 30, supports 
this critical behavior. The PF transition line ends at the 
multicritical point M located at Tm = ctai ~ 0.784 along 
the N line. For comparison, we mention that, unlike the 
gauge-glass model, the quenched disorder is relevant at 
the PF transition line of the 3D ± J Ising models, giving 
rise to a new universality class which is the same of the 
3D randomly dilute Ising systems [s^ [3^ . 

A low-T glassy phase appears for sufficiently large dis- 
order, i.e., cr > aM separated by a finite-T PG transition 
line from the paramagnetic phase. A reasonable hypoth- 
esis is that the glassy critical behavior is universal along 
the PG transition line, up to cr —)> oo (i.e., the model 
with uniform disorder distribution), and belongs to the 
3D gauge-glass universality class. This is supported by 
MC simulations for cr = 4/3 and for the uniform random- 



3 



shift distribution, up to lattice size L = 20, which pro- 
vide clear evidences of finite-T transitions in both cases 
and the universality of their glassy critical behaviors. 
Moreover, their FSS analyses give Tc = 0.475(10) and 
Tc = 0.46(1) respectively for a = 4/3 and tr oo, and 
the estimates v — 3.2(4) and 77 — —0.47(2) for the uni- 
versal exponents describing the critical overlap correla- 
tions. These critical exponents may be compared with 
those of the glassy transition in Ising-like spin glasses, 
where [1, ES SH = 2.45(15) and 77 = -0.375(10). 

The PF and PG transition lines meet at the critical 
point M along the N line, see Fig. [TJ Actually, M is 
a multicritical point, characterized by a magnetic-glassy 
multicritical behavior with two relevant perturbations (in 
the absence of external fields). The tangents at M of 
the transition lines are parallel to the T axis. Indeed, as 
proved in Ref. [l^, a^i is an upper bound for the values of 
a where the ferromagnetic phase can exist. A FSS analy- 
sis of MC simulations along the N line, up to L = 20, lo- 
cates the point M at Tm — cfm — 0.7840(2), and provides 
the estimates yi — 0.93(3) and 7/2 = 0.56(3) for the RG 
dimensions of the relevant perturbations, thus a crossover 
exponent (j) = ?/i/?/2 = 1-'''(1)- Moreover, the exponent 77 
associated with the spin and overlap correlation functions 
is 77 = —0.121(1). Again, these critical exponents may 
be compared with those at the multicritical point of the 
phase diagram of the 3D ± J Ising model where the PF 
and PG transition lines meet, which are [i^l j/i = 1.02(5), 
7/2 = 0.61(2), (f) = 1.67(10) and 77 = -0.114(3). 

A transition line separating the low-T ferromagnetic 
and glassy phases starts from the multicritical point M 
toward the T = axis, see Fig. [TJ The order parame- 
ter is provided by the magnetic variables and their cor- 
relations, which become effectively paramagnetic in the 
glassy phase. We present FSS analyses of MC simula- 
tions at fixed T < Tm up to L = 12. They show that the 
FG transition line runs almost parallel to the T axis; it 
is slightly reentrant, consistently with the fact that aM 
is an upper bound for the presence of ferromagnetism. 
For example, we find CTc = 0.777(2) at T = 0.376. The 
magnetic critical behavior varying T is compatible with 
a universal critical behavior along the FG transition line, 
with critical exponent v = 1.0(1). 

The main features of the phase diagram and the uni- 
versality classes of the magnetic and glassy critical be- 
haviors at the different transition lines are expected to 
be largely independent of the detail of the distribution. 
For example, they are expected to apply to the case of 
random shifts with Gaussian distribution. 

The paper is organized as follows. Sec. |TT] provides the 
definitions of the quantities considered in our work. In 
Sec, mil we discuss the phase diagram at low disorder and 
the critical behavior at the PF transition line. In Sec. IIVI 
we focus on the the phase diagram at large disorder, i.e., 
large values of cr, where the low-T phase is glassy, and 
study the glassy critical behavior at the PG transition 
line. In Sec. |V] we study the multicritical behavior at 
the point M of the phase diagram, where the different 



transition lines meet. Finally, in Sec. IVII we investigate 
the FG transition line, which runs from M down to the 
T = axis. 



II. NOTATIONS 

We consider the gauge-glass model ([Ij defined on cubic 
lattices of size with periodic boundary conditions. We 
define the magnetic correlation function 



G{x -y) = [{ipx -ipy)] 
and the overlap correlation function 

where is the overlap variable defined as 

using two copies Tpi^"^ and Tp3^^ for the same disorder con- 
figuration. The angular and square brackets indicate the 
thermal average and the quenched average over disorder, 
respectively. We define the magnetic and overlap suscep- 
tibilities as 



(3) 
(4) 
(5) 



(6) 



and the magnetic and overlap second-moment correlation 
lengths 



^2 ^ G(0) - G'(gn.in ) ,2 _ 
9min'-^(9min) 



' r (n ■ \ ' 

"imin^oVyrnmy 



where q^nin = (27r/L, 0, 0), q = 2sin(7/2. 

We also consider quantities that are invariant under 
RG transformations in the critical limit, such as the ra- 
tios 



and the cumulants 



(8) 



C/4 



[{WW 



and 



where ^ = V'x and = Y^xlx- Finally, given two 
replicas of the system with spins T/ji^"* and T/ji^"* , we con- 
sider the quantity fH [28j 

/„=/3^[(/W)(/(2))] (11) 
where (3 = l/T, and 

/«^l5]lm7/;«C/_+e,V'i!ie.. (12) 

X 

is the derivative of the free energy with respect to a twist 
along one direction ei of the lattice. 
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III. THE PARAMAGNETIC-FERROMAGNETIC 
TRANSITION LINE 

The random shifts of the gauge glass model ^ van- 
ish when cr — > 0, thus recovering the pure cubic-lattice 
nearest-neighbour XY model, which undergoes a contin- 
uous transition at Txy — 2.201842(5) between the high- 
T paramagnetic phase and a low-T ferromagnetic phase 
with long-range order. The critical behavior belongs to 
the 3D XY universality class [11], characterized by the 
symmetry U(l), which also describes transitions related 
to the formation of Bose-Einstein condensates in interact- 
ing quantum particle systems, the superfluid transition in 
"^He, transitions in easy-plane magnets, etc... The critical 
exponents, which determine the asymptotic behaviors of 
the critical correlations, are [13, HI] vxy = 0.6717(1), 
•q^Y = 0.0381(1), axY = 2 - Suxv = -0.0151(3), etc... 

The low-T ferromagnetic phase is expected to be stable 
with respect to the presence of a weak disorder respect- 
ing the global U(l) symmetry, such as that arising in the 
gauge glass model for nonzero values of cr, see also p^ . 
analogously to the 2D gauge glasses (where actu- 
ally we have a low-T quasi- long-range order phase), and 
the 3D ± J Ising model [39] where the disorder does not 
break the global Z2 symmetry. Therefore we expect that 
a paramagnetic-ferromagnetic (PF) transition line starts 
from the pure XY critical point P = (cr = 0,T = Txy), 
where the relevant symmetry is still U(l). 



A. Irrelevance of the disorder at the PF transition 

The critical behavior along the PF line can be inferred 
by studying the relevance of the RG perturbation induced 
by the random shifts at the pure 3D XY fixed point. We 
argue that, since the global U(l) symmetry is maintained 
in the presence of random shifts, the leading RG per- 
turbation induced by the random shifts gets effectively 
coupled to the energy density at the PF transition line, 
as in the case of randomly-dilute spin models [H, lil j. 
This implies that the relevance of the disorder is related 
to the sign of the specific-heat exponent of the pure sys- 



tem [4J, |45| . If it is positive, like the case of 3D Ising-like 
models, the disorder provides a relevant perturbation, 
which changes the asymptotic critical behavior, giving 
rise to a new universality class ,39]. On the other hand, 
if the specific-heat exponent is negative, like the case of 
the 3D XY models where axY = —0.0151(3) (but also 
for any 3D model with symmetry 0{N) and N > 2 [38|). 
the disorder is irrelevant in the RG sense, which implies 
that the critical behavior belongs to the same universality 
class of the pure system, and the asymptotic power-law 
behaviors remain unchanged. 

However, disorder gives rise to scaling corrections 
which gets suppressed very slowly and are absent in the 
pure system. Analogously to randomly dilute 3D XY 
models [s^, they are 0(^~"'*) where is the critical 
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FIG. 3: MC data of at a = 0.35. 



length scale and 



Wd = 



axY 

VXY 



= 3- 



VXY 



= 0.0225(5), 



(13) 



to be compared with the leading scaling correction of 
pure XY systems which are 0(C-") with [H HI] w = 
0.785(20). 



B. Monte Carlo simulations at a = 0.35 

We support the above scenario by a FSS analysis of MC 
simulations at cr = 0.35, up to lattice size L — 30. In the 
simulations we use both Metropolis and microcanonical 
local updatings. One single step of the Monte Carlo up- 
date was composed by one Metropolis sweep followed by 5 
microcanonical sweeps. See Ref. M for more details. We 
average over a large number of samples, from Ng — 16000 
for L = 6 decreasing to Ng — 3000 for L = 30. The equi- 
libration is carefully checked by monitoring the MC time 
evolution of the observables which we consider. 

The FSS of the RG invariant quantities, such as i?^, 
C/4, and J/l, show a clear evidence of a continuous 
transition, see, e.g., Fig.[3l A standard FSS analysis gives 
Tc = 1.7103(3) from the crossing point of their data for 
different values of L, and v — 0.68(1) from their slope at 
Tc (using data for L > I/min ^ 12). This estimate of v is 
in good agreement with the exponent vxy = 0.6717(1) of 
the 3D XY universality class. Indeed, Fig.|4]shows that a 
good collapse of the MC data of is achieved by plotting 
them versus (T - Tc)/.!/""^^ with vxy = 0.6717, thus 
supporting the universality with the 3D XY universality 
class. Analogous results are obtained for the other RG 
invariant quantities. 

Given two generic RG invariant quantities Ri and R2, 
the relation 



Ri = /fli (-R2) 



(14) 



is asymptotically universal, i.e., independent of the model 
within the given universality class, apart from scaling 
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corrections. This fact provides further stringent checks 
of universahty. As an example, Fig. [S] shows a plot of 
versus i?| from MC simulations of the gauge glass 
with a = 0.35 and of the pure XY model, i.e., u = 0. 
The data appear to approach a unique universal curve 
in the large-L limit, providing a further strong evidence 
of universality. Analogous results are obtained for other 
combinations of R^, U4 and U^. 

These numerical results provide a strong evidence that 
the critical behavior along the PF transition line belongs 
to the 3D XY universality class of pure systems. The 
above-mentioned slowly decaying scaling corrections are 
best observed in the quantity U22, cf. Eq. ([9]), which is 
trivially zero in pure systems. In Fig. |6] we show U22 
versus for several values of L. The data of U22 at 
a = 0.35 are very small, U22 < 0.03, but nonzero. U22 is 
expected to vanish for L ^ 00, but very slowly, as 
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L-^-'fu.ARi), = 0.0225(5), 



(15) 



IV. THE PARAMAGNETIC-GLASSY 
TRANSITION LINE 



where f22{R{) is a universal function apart from a triv- 
ial overall normalization. Without taking into account 
Eq. p^ . the data of U22 in Fig. |4] might be considered 
as an evidence of an unexpected scaling behavior. How- 
ever, one can easily check that they are still compatible 
with a very slow suppression in the large-L limit, such as 
Eq. (1151) . Indeed, with increasing the size from L = 10 
to L = 30, the expected variation should be just a few 
per cent, which is within the typical statistical error of 
the data in Fig. A check of the vanishing large- 

L limit of U22 would require much larger lattices and 
higher statistics. A similar situation has been met at the 
PF transition line of the 2D ± J Ising model, where the 
U22 vanishes logarithmically, as shown in Ref. |46I . 



In the case of a uniform distribution of random shifts, 
which is formally represented by the limit cr 00, 
there is already a good numerical evidence, see, e.g., 
Refs. [2^, [2^ for a finite-temperature PG transition be- 
tween the high-T paramagnetic phase and a low-T glassy 
phase. The glassy phase is expected to persist when we 
consider finite large values of cr; thus, a PG transition 
line is expected for sufficiently large values of a. More 
precisely, it is expected to run from the large disorder 
limit down to the multicritical point M, see Fig. [1] at 
Ta[ = aM ~ 0.784, see the next section. Our working 
hypothesis is that the critical behavior along this PG line 
is universal, belonging to the 3D gauge-glass universality 
class. 
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A. Monte Carlo simulations 

In order to investigate the glassy critical behavior at 
large disorder, we present FSS analyses of MC simula- 
tions for (7 = 4/3 and for uniformly distributed phase 
shifts, formally corresponding to cr — ?■ oo, up to lattice 
sizes L ~ 16 and L = 20 respectively. We perform av- 
erages over lO"* disorder configurations for all the lattice 
sizes and temperatures considered. 

In these MC simulations we supplement the local MC 
updating method, used at the PF transition line, with 
the random-exchange or parallel tempering method, see, 
e.g, Ref. lU which allows us to reliably simulate small 
values of T, in particular below the N line T — a. In 
the parallel-tempering simulations we consider Nt sys- 
tems at the same value of a and at Nt different temper- 
atures in the range Tmax ^ Ti > Tmin, with Tmax ~ 0.87 
and Tmin ~ 0.37. The value Tmax is chosen so that the 
thermalization at Tmax is sufficiently fast, while the in- 
termediate values Ti are chosen so that the acceptance 
probability of the temperature exchange is at least 5%. 
Moreover, we require that one of the Ti is along the N 
line, i.e., Ti = a, where the known exact results allow 
us to check the MC code and the thermalization. The 
thermalization is further checked by verifying that the 
averages of the observables remain stable for all Ti after 
a sufficiently large number of MC steps for each disor- 
der realization. The overlap correlations and the cor- 
responding Xo and SLie measured by performing two 
independent runs for each disorder sample. In the case 
of observables requiring the computation of the disorder 
average of products of thermal expectations, such as the 
case of U^, we use unbiased estimators as explained in 
Refs. 0, S (a naive application of the disorder average 
would introduce a bias, thus a systematic error). 



B. Finite-size scaling analysis 

In Fig. [7] we show the MC data of i?| and lo, cf. 
Eq. (fTTj) . Their crossing points provide a strong evidence 
of a finite-T transition for both cr = 4/3 and a = oo. A 
large exponent v is already suggested by the slow rising 
of the slopes at the crossing point with increasing L. 

To beghin with, we address the universality issue, i.e., 
whether the transitions at u = 4/3 and cr = oo have 
the same universal critical behavior. For this purpose, 
in Fig. [5] we plot data of Jo, and C/22 versus i?| at 
CT = 4/3 and cr = 00. They appear to converge toward 
the same universal large-i limit, providing a strong evi- 
dence of universality. These results support our working 
hypothesis that the glassy critical behavior is univeral 
alsong the PC transition line from large disorder to the 
multicritical point M. 

The FSS behavior of the RG invariant quantities allows 
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FIG. 7: i?| for a = 4/3 (above), and (middle) and h 
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us to estimate Tc and v: 

R{T,L) = ,f{utL^'^) = 

^R* + ci{T- T^)L^/'' + C2{T- TcfL^''' ■ 



(16) 



where ut is the temperature scaling field, ut ^ T — T^ 
R indicates the generic RG invariant quantity, such as 
and lo, and we neglect scaling-correction terms. We 
check the stability of the fits by increasing the minimum 
size Lmin of the data used in the fit, varying the range 
of values of T around Tc fwe use self-consistent windows 
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around Tc limiting the value of (T — Tc)L^^'^ , which cor- 
responds to limiting the range of values of any R around 
R*], and the number of terms in Eq. We use the 

comparison of the results using different quantities as a 
check of the relevance of the neglected scaling corrections. 
Some results of fits of the data of R^ and 1^,, using the 
Ansatz (fT6|) . are reported in Table HI 

In the case of the uniform distribution we have data 
up to L = 20. The fits of the data of i?| are reasonably 



TABLE I: Results of fits to R* + Ci (T -T^yu'" of the 

data of and lo, with respect to variations of the range of 
T (taking data for — _R|*| < A_R| to have a self-consistent 
scaling cut, with ~ 0.7), the number n of terms in the fit 
Ansatz, and Lmin which is the minimum size allowed by the 
data. 



R ARs^ 


n 




Tc 


1/u 


R* 


xVdof 


00 i?| 0.07 


1 


8 


0.467(2) 


0.27(3) 


0.684(5) 


3.1 






10 


0.477(3) 


0.29(4) 


0.665(7) 


2.9 






12 


0.458(4) 


0.31(7) 


0.703(9) 


1.8 


0.1 


1 


8 


0.465(2) 


0.27(2) 


0.689(4) 


3.6 






10 


0.475(2) 


0.29(2) 


0.669(5) 


3.5 






12 


0.455(4) 


0.29(4) 


0.712(9) 


2.1 






14 


0.461(5) 


0.33(6) 


0.698(11) 


1.8 


0.1 


2 


8 


0.466(2) 


0.29(8) 


0.687(3) 


3.2 


0.2 


2 


6 


0.463(1) 


0.31(2) 


0.690(2) 


3.3 






8 


0.460(1) 


0.32(2) 


0.698(2) 


3.0 






10 


0.469(2) 


0.32(2) 


0.680(3) 


2.8 






12 


0.450(3) 


0.29(3) 


0.721(6) 


1.6 






14 


0.456(4) 


0.31(4) 


0.708(9) 


1.3 


lo 0.1 


1 


8 


0.467(5) 


0.26(4) 


1.74(5) 


1.4 






10 


0.472(6) 


0.27(5) 


1.69(5) 


1.5 






12 


0.488(7) 


0.33(7) 


1.53(7) 


0.8 


0.3 


2 


8 


0.463(2) 


0.38(12) 


1.77(2) 


1.6 


4/3 R^ 0.1 


1 


6 


0.486(2) 


0.31(2) 


0.669(4) 


2.7 






8 


0.477(2) 


0.31(2) 


0.685(5) 


1.5 






10 


0.478(5) 


0.28(5) 


0.683(10) 


1.8 


0.1 


2 


8 


0.476(3) 


0.31(11) 


0.686(4) 


1.0 


0.2 


2 


6 


0.479(1) 


0.33(2) 


0.682(2) 


3.1 






8 


0.471(2) 


0.33(2) 


0.695(3) 


1.4 






10 


0.471(3) 


0.31(3) 


0.695(6) 


1.6 



stable. We note some oscillations in the estimates of Tc, 
which may reflect the difficulty to estimate when the 
critical exponent v is large. From these results one may 
get the estimates = 0.46(1), l/v = 0.31(4) and = 
0.70(2), where errors take also into account the stability 
with respect to changes of L^un, from Lmin = 6 to Lmin = 
12, and the range of T around the crossing point. The 
results using lo are definitely consistent, although they 
appear less precise (also because they correspond to less 
statistics, we started later collecting data for /q). They 
suggest the estimates = 0.47(2), l/iy = 0.30(5) and 
/* = 1.7(1). The fits of the data for the cumulants f/| 
and U22 do not provide sufficiently stable results, but 
they scale consistently as shown by Fig. [HI 

In the case of the distribution with a = 4/3, whose 
data are up to L = 16, the results appear more stable. 
We obtain = 0.475(10), l/v = 0.32(4) and = 
0.69(1) from the fits of the data of R°. For a = 4/3 we 
have only few data for lo, see Fig. [SI which do not allow 
us to obtain an independent estimate of v. 

The above results for v are consistent, thus, in agree- 
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0{h?). Thus, neglecting nonanalytic scaling corrections, 
we fit the data of Xo to 

inxo -aini + &o + &i(r-r,) + ... + ci(r-r,)ii/'' + ... 

(19) 

fixing the values of Tc and v as obtained above. We 
obtain 



f] = -0.47(2) 



(20) 



for the uniform distribution and 7] = —0.46(2) for the 
distribution with a — 4/3, where the error takes also 
into account the uncertainty on and v. 

Finally, we compare our results with earlier numeri- 
cal works for the gauge-glass model with uniformly dis- 
tributed random shifts. Earlier estimates of Tc, i.e., 
Tc = 0.47(3) from Ref. [H and = 0.48(2) from Ref. HO, 
are consistent with ours. Our result v — 3.2(4) is sig- 
nificantly larger than the estimates obtained by earlier 
FSS analyses, u = 1.3(4) from Ref.[li, v = 1.39(20) from 
Ref.^, and v = 1.62(20) Ref. [10. These smaller values 
of V appear excluded by our FSS analyses. This discrep- 
ancy may be explained by the small lattice sizes, up to 
L — 12, of their data. We also mention that our esti- 
mate is quite larger than the experimental estimate at 
the vortex-glass transition in the (K,Ba)Bi03 cubic su- 
percondutor reported in Ref. [H, i.e., v = 1-0(2), which 
has been often compared with the glassy transition of 
the gauge glass with uniform distribution. Finally, our 
estimate (^01 of 77 agrees with the result rj — —0.47(7) 
reported in Ref. [2^. 



FIG. 9: lo (above) and (bottom) versus (T - Tc)L^''" for 
the uniform distribution, with Tc = 0.46 and 1/v = 0.31. 



THE MULTICRITICAL POINT ALONG THE 
N LINE 



ment with universality. We consider 

1/i^ = 0.31(4), 1/ = 3.2(4), (17) 

as our final estimate. Fig.[S]shows the collapse of the data 
of i?| and lo for the gauge glass with uniform distribu- 
tion, when they are plotted versus (T - Tc)L^/'', with 
Tc — 0.46 and 1/j/ = 0.31. Note that is quite large, 
but still comparable with the value v = 2.45(l5)of the 
glassy transition in 3D Ising-like spin glasses [1, |4l| . 
We should also say that the accuracy and the relatively 
small lattice size of the available data do not allow us 
a robust control of the neglected scaling-corrections, as 
it was achieved for the Ising-like spin glass models [49| . 
Therefore, further numerical work is required to substan- 
tiate the above results. 

We compute the exponent 77 associated with the over- 
lap correlation ^ , by analyzing the data of the overlap 
susceptibility Xo, which is expected to behave as 0| 

Xo = uW-'^"fiutL'^'') . (18) 

Here u/j is related to the external overlap scaling field 
associated with the overlap variables by = huh{T) + 



In the case of the cosine distribution ([2]), the N line 
plays an important role in the phase diagram, because 
it marks the crossover between the magnetic-dominated 
region and the disorder-dominated one. We conjecture 
that the PF and PG transition lines meet at a multicriti- 
cal point M, which coincides with the critical point along 
the N line, analogousl y to the phase diagram of the 3D 
±J Ising model 0, [13,113 ■ This is also suggested by the 
fact that the location of the critical point along the N 
line provides a bound on the disorder parameter a where 
a ferromagnetic phase can exist [lB|- We also mention 
that the critical point along the N line shows a multi- 
critical behavior also in the 2D ± J Ising model and 2D 
gauge glass, see, e.g., Refs. [l^ [sol. [sTI. even though these 
models do not have a low-T glassy phase, thus a finite-T 
PG transition line. We shall support the above scenario 
by a FSS of MC simulations of the 3D gauge glass along 
the N hue. 



A. The Nishimori line 

The cosine distribution ^ lends itself to some exact 
calculations along the N line [H, [i3|- For example, the 
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energy density is^E = -3/i(l/r)/Jn(l/T) when T = 
a, where are modified Bessel functions. Moreover, the 
spin and overlap correlation functions are equal: 

[{i'.i'y)] = m.i^y)n (21) 

Along the N line we also have = i?| and U4 — U^. 

As proved in Ref. [IB], the critical value (tm of a along 
the N line is an upper bound for the values of a where 
the ferromagnetic long-range order can exist. Therefore, 
at the critical point M = {aM,TM) the tangent to the 
transition line limiting the ferromagnetic phase must be 
parallel to the T axis; moreover, the critical value ao 
where the ferromagnetism disappears at T — must sat- 
isfy ao < (JM- 

The N line T = cr is also characterized by an exten- 
sion of the replica symmetry. The quenched average over 
disorder, which implies the average of the free energy, 
can be formally reconstructed by introducing n replicas 
of the system, taking the limit n ^ after the disorder 
average. We write the disorder average over the partition 
function of n replicas as 




j ni^V'^")] exp[(l/T) ^ cos(e) - - A.,y)] 

(xy) 

Let us perform the gauge transformation 

i^x ^ e'^^ , Uxy ^ e-'^^ Uxye"^y . (23) 
We obtain 




I llm^-^] exp[(l/T) ^ cos(e^ - e^"^ - A^y)]. (24) 

Ji (xy) 

Since a further integration with respect to (px = e**^^ gives 
only rise to a trivial constant factor, we have that the 
gauge variables (j)x correspond to another identical replica 
when T = a, extending the original replica symmetry. 

B. Scaling behavior at the multicritical point 

As a working hypothesis, we assume that the critical 
point along the N line is a multicritical point, analogously 
to the case of the 3D ± J Ising model [4^ . We derive some 
predictions which are then verified by a FSS analysis of 
numerical MC simulations. 

In the absence of external fields, the critical behav- 
ior at the multicritical point M is characterized by two 
relevant RG operators. The singular part of the free en- 
ergy averaged over disorder in a volume of size L can be 
written as 

F,i^g{T,a,L)=L-''f{uiLy\u2Ly\{u,Ly^}), (25) 



with i > 3, where yi > 2/2 > 0, < for i > 3, Ui 
are the corresponding scaling fields, and ui = M2 = at 
M. In the infinite- volume limit and neglecting subleading 
corrections, we have 

i^.ing(r,a) = |7.2|''/^V±(«i|«2|-^) (26) 

around M, with (jj — yi/y2 > 1, where the functions 
f±{x) apply to the parameter regions in which ±U2 > 0. 
Close to M, all transition lines correspond to constant 
values of the product ui|m2|~'^ and thus, since > 1, 
they are tangent to the line ui — 0. 

The relevant scaling field ui and U2 can be inferred 
by using the following facts: (i) since gm is an upper 
bound for the values of a where the ferromagnetic phase 
can exist, the transition line at M must be parallel to 
the T axis; (ii) the condition T = cr at the N line is RG 
invariant, because it is protected by the extension of the 
replica symmetry, as shown in Sec. IV Al We therefore 
have 

Ul=(7 - (7M + ■■; (27) 

where the dots indicate nonlinear corrections, which are 
quadratic in Act = u — aM and AT = T — Tm, so that the 
line Ml = runs parallel to the T axis at M. Moreover, 
we choose 

U2 = T-a, (28) 

so that the N line corresponds to U2 = 0. 

These results give rise to the following predictions for 
the FSS behavior around M. Let us consider a RG invari- 
ant quantity i?, such as R^, C/4, C/22, defined in Sec. [Til 
In general, in the FSS limit R obeys the scaling law 

R = n{uiLy\u2Ly\{uiLy^}), i>3. (29) 

Neglecting the scaling corrections which vanish in the 
limit L — >■ 00, we expect 

R ^ R* + bnUiLy^ + b2iU2Ly^ + . . . . (30) 

Along the N line, the scaling field U2 vanishes, so that we 
can write 

Rn = r* + hiuiiy' + h2iulL^y' + (31) 

where the subscript N indicates that R is restricted to the 
N fine. Differentiating Eq. with respect to ,3 = 1/T, 
we obtain 

R' = biiu^LV^ + b2iu'^Ly^ + • • • (32) 

If Eq. (gZl) holds, then u\ = 0(T - Tm), so that the 
leading behavior along the N line is 

R'n = ^21^2-^^^ + • • • (33) 

The magnetic susceptibility along the N line behaves as 

XN = eL^-" (1 + eiu^Ly^ +•••)• (34) 

Note that there is only one "q exponent which charac- 
terizes the critical behavior of both the magnetic and 
overlap correlation functions, since they are equal along 
the N line. 
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FIG. 10: i?5 along the N line T = a. 





C. MC results 

In the following we present a FSS analysis of MC sim- 
ulations along the N line. The MC algorithm was a mix- 
ture of standard Metropolis and microcanonical updates, 
as in the MC simulations at the PF transition line re- 
ported in Sec. IIIII In order to locate the multicritical 
point we simulated several temperatures ranging from 
T ^ 0.773 to T = 0.797, for lattice sizes 6 < L < 20. The 
number of disorder configurations ranged from 4 • 10^ to 
2 • 10^ for the largest lattices. This large statistics was 
necessary to achieve a convincing evidence of the multi- 
critical nature of M. 

The MC data of are shown in Fig. \M There is 
clearly a crossing point at T w 0.784. Analogous results 
are obtained from U4 and U22- In order to estimate Ta/ 
and j/i, we fit the RG invariant R to 

R = R*+ai{T- TM)Ly' + a2{T ^ TM fL^"' + .... (35) 

Note that this functional form relies on the property that 
U2 = along the N line. Otherwise, an additional term 
of the form (T - TM)Ly^ should be added. We also ne- 
glect scaling corrections which are 0{Ly^) with 7/3 < 0. 
We check their relevance by comparing results from the 
analyses of different quantities. 
Our best estimates are 

Tm = <TAf - 0.7840(2), (36) 
2/1 = 0.93(3), (37) 

and also i?| = 0.5594(4), Ul = 1.226(1) and U^^ = 
0.128(4). The errors take into account the stability of the 
results with respect to the minimum size Lmin allowed in 
the fits (typically from Lmin = 6 to imin — 12), the range 
of values of T around [we use again self-consistent win- 
dows around Tc limiting the values of (T — Tm)^^^], the 
number of terms in Eq. (j35p . Scaling corrections are ap- 
parently small. We find some evidence of scaling correc- 
tions only in the analysis of C/22, but they appear to decay 



quite fast, suggesting a relatively large scaling correction 
exponent, i.e., 2/3 w —2. 

The derivative R' with respect to (3 can be estimated by 
computing appropriate connected correlations of R and 
the Hamiltonian in the MC simulations along the N line. 
According to the multicritical scenario outlined above, 
R' is expected to behave as L^^ at Tm with 1/2 < Ui- In 
order to determine 7/2, we fit the data to 

\nR' ^a + y2\riL + h{T -TM)Ly\ (38) 

keeping Tm = 0.7840 and yi = 0.93 fixed. We obtain 

2/2 = 0.56(3) (39) 

from R'^ (the error includes the uncertainty on Tm and 
2/1), while the data of C/4 are not sufficiently precise to 
provide a stable result. The above estimate of 2/2 nicely 
support the multicritical scenario, since it shows that 
2/2 < 2/1. Therefore the crossover exponent, cf. Eq. (|26p . 
is 

0=^ = 1.7(1). (40) 

2/2 

We determine the exponent 77 from the FSS of the ratio 
Z = x/^^ i"'' at the critical point. We fit its MC data 
along the N line to 

In Z = a - r; In L + b{T - TM)Ly^ + c{T - Tm) (41) 

where the last term takes into account possible analytic 
terms coming from the scaling fields [4], analogously to 
Eq. (dH]). We obtain 

= -0.121(1) (42) 

It is worth noting that the above estimates of the mul- 
ticritical exponents are quite close to those found for the 
3D ± J Isingmodel, at the multicritical point along its N 
line, where [13 2/1 = 1.02(5), 2/2 = 0.61(2), = 1.67(10) 
and 77 = -0.114(3). 



VI. THE FERROMAGNETIC-GLASSY 
TRANSITION LINE 

The EG transition line runs from the multicritical point 
down to T 0, with Uc < (Tmj since ctm provides a bound 
for the region where the magnetic order can exist. The 
order parameter along this transition is provided by the 
magnetic variables and their correlations. 

Ref. 16 argues that this transition line runs parallel to 
the T axis, thus CTc = cm for T < Tm- An analogous 
prediction for the 2D ±J Ising model turns out to fail, 
although it provides a good approximation, because the 
low-T transition line where ferromagnetism disappears 
turns out to be almost parallel and only slightly reen- 
trant, see, e.g., Refs. HH, [H and references therein. 

We investigate this issue by numerical MC simu- 
lations up to lattice sizes L = 12, using the same 
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FIG. 11: i?e (below) and (above) at T = 0.376 across 
the FG transition line. The vertical dashed lines show the 
estimate of Tm, i.e., Tjvj — 0.7480(2), with its uncertainty. 



over a large number of of disorder configurations: = 
16000, 8000, 6000, 4000 respectively for L = 6, 8, 10, 12. 
Fig. [n] shows the data of and i?| at T = 0.376. The 
set of data of for different lattice sizes show a crossing 
point, confirming the existence of FG transition. On the 
contrary, the data of i?| do not show crossings, which 
may reflect the fact that such transition separates two 
ordered phases with respect to the overlap variables. 

Note that the crossing points of the i?j data appear 
to cluster at a value of a which is slightly smaller than 
aM = 0.7840(2). Indeed, fitting them to 

R^=Rl + ai{a-ac)L^/'' + a2{a-acfL^/'' + ..., (43) 
we obtain 



cTc 0.777(2), = 1.0(1), 



(44) 



and i?| = 0.73(1). These estimates should be taken with 
some caution, in particular that of due to the small 
size of the available lattices, which does not allow us to 
perform stringent checks of stability. The analysis of the 
data of U4 gives consistent results, but less precise. We 
also mention that an analogous FSS analysis of the data 
at a larger temperature T = 0.437, but still smaller than 
Tm, gives ct^ = 0.782(2), ly = 1.1(1) and i?| = 0.71(1), 
which support the universality of the magnetic critical 
behavior along the FG transition line. 

According to the above results, the critical values of 
the disorder parameter at T < Tm are very close but 
smaller than cm = 0.7840(2), indicating a slight reen- 
trant transition line. 
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MC method employed at the glassy transitions. Us- 
ing random-exchange MC simulations, we collect data 
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